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EFEITO DA VARIAÇÃO AXIAL DA TEMPERATURA OU DO 
FLUXO TÉRMICO NO COEFICIENTE DE TRANSMISSÃO DO 
CALOR POR CONVECÇÃO FORÇADA NUM TUBO 


RESUMO 


Discute-se a importância das condições na superficie 
de permuta — temperatura ou fluxo térmico impostos — no 
valor do coeficiente de transmissão do calor por convecção, 
num tubo, os quais podem diferir substancialmente, Utili- 
zando propriedades das equações fundamentais mostra-se 
como se podem calcular os coeficientes quando as condições 
impostas variam axialmente desde que seja conhecido o seu 
valor (por via teórica ou experimental) para condições 
uniformes, 


INTRODUÇÃO 


Em todos os textos sobre transmissão de calor, 
com algumas preocupações de rigor, se faz a 
distinção entre coeficiente de transmissão por 
convecção forçada (h) quando a temperatura na 
superfície que delimita o escoamento é uniforme 
do caso em que é uniforme o fluxo térmico. Sob 
o aspecto matemático, o primeiro caso corres- 
ponde a impor como condição ao limite uma 
temperatura conhecida; a segunda a impor como 
condição o valor da sua derivada. Estas formas 
de condições dos limites escolhem-se pela sua 
simplicidade e também porque na prática corres- 
pondem muitas vezes a dois limites de funciona- 
mento: um permutador em que um dos fluidos 
muda de fase aproxima-se, relativamente ao outro 
fluido, do caso de temperatura na superfície cons- 
tante; num permutador em contracorrente em 
que os fluidos têm igual G.C,, realiza-se apro- 
ximadamente o caso dum fluxo térmico constante 
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por J.J. DELGADO DOMINGOS 


Prof. do L.S.T. 
Núcleo de Estudos de Engenharia 
Mecânica (do 1,A.C,) 


SYNOPSIS 


The importance of boundary conditions — uniform 
temperature or heat fluz —on the heat transfer coeficient 
by forced convection in a tube is discussed, 

Starting from the properties of the fundamental equations 


for heat transfer it is shown how the heat transfer 


coeficient can be found for axially varying boundary 
conditions when its value is known (experimentally or 
analytically) for unijorm temperature or heat Aux, 


entre eles porque a sua diferença de temperatura 
também permanece constante. Existem porém 
muitos casos intermédios em que o projectista 
apenas se baseia na sua intuição e se escuda em 
coeficientes de segurança de que resultam quase 
sempre permutadores menos económicos, porque 
calculados com valores demasiado baixos para o 
coeficiente de transmissão. Noutros casos, como 
por exemplo em Reactores Nucleares, correm-se 
riscos que envolvem a segurança do núcleo por- 
que há deficiente previsão não apenas da tempe- 
ratura máxima como da posição em que a mesma 
se verifica nos canais de arrefecimento. 


0.1. O exemplo mais simples (mas nem por 
isso dos menos importantes) que se pode consi- 
derar é o da transmissão do calor num tubo cir- 
cular, em que o escoamento é uniforme (*). 


(*) ou regime estabelecido: o perfil de velocidade não 
varia axialmente, 
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Se além do escoamento ser uniforme e laminar 
o regime térmico se encontra estabelecido e as 
propriedades físicas são constantes o número de 
Nusselt assume os valores: 


N —hd 3,658 Se a temperatura da 
k parede é uniforme 


hd Se é uniforme o fluxo 
térmico imposto na pa- 
rede 


Como se verifica o coeficiente de transmissão 
no primeiro caso é 16/ inferior ao do segundo, 
o que, mesmo em muitos casos de prática corrente 
se não pode considerar desprezável. 


0.2. Os valores anteriores correspondem ao 
regime térmico estabelecido. Como é sabido, o 
coeficiente de transmissão é máximo na secção 
em que se inicia o aquecimento, tendendo assimp- 
tóticamente para os valores que correspondem 
ao regime estabelecido com o deslocamento para 
juzante; o mesmo sucede com a diferença entre 
os N, para os casos considerados. 

Como é óbvio, sob o aspecto aplicado é impor- 
tante saber calcular o N, quando nem a tempe- 
ratura nem o fluxo imposto na superfície de 
permuta são axialmente uniformes. Basta, aliás, 
obter a solução para um dos casos pois uma 
escolha conveniente da variação de valor imposto 
permite realizar não só a condição ao limite 
correspondente ao outro caso, como a todas as 
situações de interesse. 

Essa generalização é o objectivo do presente 
trabalho. Os resultados aplicam-se tanto ao 
regime laminar como ao regime turbulento e, 
embora utilizem propriedades das equações fun- 
damentais da mecânica dos fluidos e transmissão 
do calor não exigem o conhecimento da solução 
analítica para o perfil de velocidades e tempera- 
tura. Basta o conhecimento (experimental, even- 
tualmente) do número de Nusselt (*) em função 
da distância à secção onde se inicia o aqueci- 


mento. 


1 MODELO FÍSICO 
1.0. Admitiremos que: 


a) O escoamento é uniforme pelo menos a 
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partir da secção onde se inicia a transmis- 
são de calor. 

b) As propriedades do fluido não variam 
axialmente (embora possam não ser unifor- 
mes na secção transversal). 

c) A temperatura do fluido é uniforme antes 
de se iniciar a permuta de calor. 

d) A condução axial no fluido (**) é despre- 
zável, tal como na superfície de permuta. 
Não se tem em conta a dissipação viscosa. 

e) O regime é estacionário. 

f) Existe simetria axial quanto às condições 
impostas na superfície. 

Fazendo coincidir o eixo dos xx com o eixo do 
escoamento e designando por x(y,z) a velocidade 
do fluido, a equação de conservação da energia 
escreve-se [1] atendendo às hipóteses anteriores. 


() (-) (-) 
1) ee Equação Ada dica E ) 
dx dy dy dz dz 


Sendo O =T-T.. 


Para o tubo de secção circular, utilizando coor- 
denadas cilíndricas a equação anterior escreve-se: 


Devendo a solução de 1) ou 1') obedecer às 
condições dos limites: 


2) O=o parax=o 


) 6) 
3) k a =], Sobre a superfície se o fluxo 
dn 


é imposto (***). 
ou 


4) O =€), Sobre a superfície de permuta se a 
temperatura é imposta(*). 

As condições 3) e 4) correspondem os valores 
assimptóticos referidos em 0.1 para o caso do 


(*) ou de h, como é óbvio. 
(**) O que é admissível para P, > 100. 
(***) Para r=R no caso do tubo 
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tubo em regime laminar estabelecido ou os 
conhecidos para regime turbulento [2/(*). 

Para o regime térmico não estabelecido admita- 
mos que é conhecido (p. ex. experimentalmente) 
o número de Nusselt em função de x: 


Nuxq = f(x) para o caso do fluxo imposto na 
superfície (condição 2) e 3)) mas 
constante (**) 

Nut = f, (x) para o caso de temperaturas impos- 
tas mas constantes (condição 2) e 4) 


O problema consiste em saber o valor de N,. 
quando, na superfície, q ou T variarem axialmente, 


dO 

» Emi) xa 
dn 

4º) e =qd(x) x:B 


e que designaremos respectivamente por Nega e 
NuxT E) E 


2. SOLUÇÃO DO PROBLEMA 


2.1 Nas hipóteses do modelo, 1) (ou 1”) é 
linear e homogênea. 


Além disso, sendo de primeira ordem em — 
e O< x < co demonstra-se facilmente [1,3] que: 
se O (x,y) é a solução correspondente a 1,2 e 3) 


a solução correspondente a 1) 2) e 3') é dada 
por (1) 


E mm dO 
t= — — mm dz 
di q). 9% ai E 


ou 


2 pão) 
t= — —— d 
) = Ff q (2) bs; T 


a qual é também, «mutatis mutandis» a corres- 
pondente a 4”). 
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2.2. Por definição 


(E) 
on /B q k 
7 h == = =N— e ea 
Rd Go 4º | 
sendo 

[rue Ods 
8) On = = — [que Ods 

ouc, ds º 


a temperatura «adiabática» ou de «mistura» na 
secção do escoamento [3]. 

Quando o fluxo na superfície varia segundo 
3") será 


a) h* =— 2.969) + — N* k 
xq d 


a o. em 


m B 


mas, por balanço térmico 
10) O (o) = [" LO qu ("açao 
o Cc 5 o 


e O*, é dada por 5) ou 6) pondo P=B 


11) 4%, = DL IN otica (SN de 
]q Jo dx (z,P) 


Por hipótese, N... é conhecida pelo que, da defi- 


Uxg 


nição 7) 

d 
12) 0, = 0 — 7 

k N 

Uxq 
e como 
13) O, = [ A gi 4 

o E * 


resulta substituindo 13) em 12), 12) em 11) e 11) 
e 10) em 9) 


*) Note-se que os valores assimptóticos são indepen- 
dentes da condição imposta em x = O. 

**; Usamos o índice x para relembrar a variação axial 
de N, e os índices q ou T consoante se trata do 
fluxo ou da temperatura imposta. 

***) Supõe-se q (x) ou 7; (x) continua em x. Se o for 
apenas por troços resultam termos aditivos em 5) 
ou 6) que se determinam imediatamente atendendo 
às propriedades da transformada de Laplace [4). 

“***) Como é óbvio, h,, corresponde a impor o valor 
do numerador, h, a impor Op. 
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14) hºq = Nºua d as EY, 1) dz 
|. ci a [um 
e 
, o 
15) DE aura Ee 


x d 1 dz 
f D(x — Z) E q ) se 


o d uxa 


2.3. Considere-se agora o caso da tempera- 
tura na superfície ser variável, isto é, a determi- 
nação de N*, ., quando N, + é conhecido. 


Da definição de h;r (ou N,.+) resulta que 


o h pdx 
16) 9 ii Om = 9, i É 


e que 
1 x 
Otam(— Ordem qt dt=— 
CS o 
=) hat" jdz 
É » Eai m 


do, 
dx = to )=g"ar 


17) 
o pa d 6) 
O*(x,P)= m(x—T [52] ZP 
a), trt Dz) (2P) 
pelo que 


1 - 
a E” aaa 
i + 1 x q E) 
E us —T)(—— | dz ds; 
E Ja Ff: pr (x Hm) : 


ad x 
permutando J e f , € atendendo à defini- 


/ o 


ção de O. 


18) PA (xj= E. IN Pp (x—Z) e) dz 
bp ” y 


donde, por 16) 
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19) 0%, (x) = el. her ?r (x — 2) dz 


e atendendo à definição de h* .. 


d x 
die h,;r?y(x—t) dz 
RS do DE 
J hor? (x—Z) dx— Cor 
e 
d Px 
e! Nar ?y(x— 2) dz 
21) Ni iii = a Eid 


f Naa Pr(x—T) dis Co, 


2.4 — O problema enunciado encontra-se pois 
formalmente resolvido pelas expressões 14), 15), 
20), e 21). Na sua solução admitiu-se que Non 
e N,xr eram conhecidos por via experimental ou 
teórica. 


Deixando para discussão ulterior as várias 
correlações teóricas até hoje propostas, afigu- 
ra-se de interesse mostrar como as mesmas 
podem ser obtidas teóricamente, o que será 
objecto do capítulo seguinte. 


3. SOLUÇÕES EXACTAS PARA N,. EN 
EM REGIME LAMINAR NO CASO DUM 
TUBO DE SECÇÃO CIRCULAR [5,6] 


3.0 O problema consiste na solução de 1”) 
sujeito a 2) e a 3) ou 4) consoante se tratar de 


IN st ou NuxTº 


Considerem-se as variáveis adimensionais 


T,—T 
O! = 
LT, 
= 
R 
22) 
u 
* =— 
u U. 
x 
x* = E 
R, A 


TECNICA Nº 385 


Introduzidas em 1”) resulta: 


9 O+ 
= (RP) ox 


e Tr 


dt* ut  doO* l 


o + 


] 
dr*” 7” o 2 dx* 
sendo ut=2(I-—r*) a solução da equação do 


movimento para escoamento uniforme, a equação 
anterior pode escrever-se (*): 


d“O* 1 d0* JO 
+— — = (1-1) — 


dr*? rt dr* dx* 


23) 


e 2) 3) e 4) transformam-se respectivamente em 
24) O0*==1] para x*=O 


dai Ri 


25) ( 1) UT —T) 


para 1* =| 


A solução do problema pode obter-se pelo 
método de separação de variáveis [7], admitindo 
que 
27) 0% (rt, x*)=W(r* ).X(x*) 
pelo qué y e X devem satisfazer a 


28) X'+1,X=0 


ara À 


rf 


AX d (I-r*")=0 


A solução de 28) é 


Sendo a, uma constante arbitrária e ), o valor 
próprio determinado por 29) e 25) ou 26) con- 
soante a condição ao limite adoptada. 

Pode demonstrar-se que os 7, são em número 
infinito e que a cada ), corresponde > solução 
de 29), os quais constituem um sistema completo. 

Deste modo, a solução final é da forma 


O (xt,rt)= 3 abilrt)e - 


= O 


(*) Despreza se o último termo se R, P, > 100, pelas 
razões já apresentadas. 
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Sendo os a, os coeficientes do desenvolvimento 
de O* (r*,o)==1 em termos dos 4. 


3.1. Designando por G, o valor assumido por 
ba Tt=] 


G, = Ya (1) 


no caso em que 26) é a condição no limite, 
obtém-se que 


31) ok “o E 1 dá 
q— 4 2 Ge Es) 
oo seção x* 
32) 0! =8 E — e 
E 
n 
2 
E — à x* 
2G, 
33) Naxt o. | 
o G —A x* 
25 O 


os valores de ), e G, calculadas por Sellars et al 
[5] são os seguintes; 


| 
n | hn? Ga 
0 7312 0,749 
1 44,62 0,544 
2 113,8 0,463 
3 215,2 U,414 
4 348,5 0,882 (*) 


Para x! >> 0,1, apenas o primeiro termo dos 
desenvolvimentos é importante e obtém-se que (**) 


) N , e 
34 - 
uT 
2 


= 3,658 


que corresponde ao regime térmico estabelecido. 


3.2. Anaálogamente, obter-se-á para o caso de 
fluxo imposto, condições 


[9.6] 
35) N a E 
44 2 1 


Sendo 


H,= (1) em que satisfaz neste caso a 25). 
Os valores de à, e H, calculados por Sparrow 
et al [6] são 


m Nai? Hm 

| 25,68 7,880 x< 10º 

9 | 83,86 2.058 >< 107º 

3 174,2 0,901 >< 10º 
296,5 0,487 >< 10 

5 450,9 0,297 xt 10º (*) 


=, 52 


4 
(*) Param >5 1,=+4m E | Hm = 0,358 Am 


8 —1 /3 
(*) Paran>2, 1 =4n+ E ; Gy = 1,01276 +, 


(3 x! >> 0,1 define aproximadamente a região de 
entrada que será dada por x= 0,05R, P,.d 


4 — EXEMPLO 


Considere-se o caso do canal de arrefecimento 
dum reactor nuclear, onde o fluxo térmiço varia 
sinusoidalmente na direcção axial segundo: 


dp 11% ) 
— — == sen 
Jpmax L 
sendo L o comprimento do canal. 
Usando a expressão exacta para N, «q dada por 
35) e introduzindo em N*,., dado por 21) obtém- 


-se para a diferença de temperatura entre a super- 
fície e o fluido de arrefecimento [2] 


Tau E 4 
GT, To | R (1— cos (Ax*)) r 
] 
+ sen (A x!) 9 —— — e 
O Ha (O, + A) 
] 
— À Cos (A x*) Dc rm 
2 2 A (1º+ A?) 
a 


e 


Gi 0 q 4 ) 
AMA +A) 


mam 


A seguir representa-se de modo aproximado a 
forma de solução, tal como a temperatura da 
superfície fazendo o cálculo para um valor cons- 
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tante de N,, , correspondente ao regime estabele- 
cido. 


Como se verifica, o cálculo habitual com 
N, q = constante subestima a temperatura máxima 
da superfície e a posição onde a mesma se veri- 
fica. Obviamente, o valor do erro dependerá de 
Jmaxr mas o exemplo mostra como pode ser 
importante a consideração dum cálculo mais 
exacto, pelo menos em condições extremas de 
funcionamento. 
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SIMBOLOGIA 


B fronteira do domínio — intersecção da 
superfície de permuta por um plano 
normal ao escoamento 
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P — ponto genérico em (x,y) 
d — diâmetro hidráulico 

R — raio exterior do tubo 

r — coordenada radial 

n — normal exterior a B 


h — coeficiente de transmissão de calor por 
convecção 


k — condutibilidade térmica do fluido 


k, — condutibilidade térmica total (molecular + 
+ turbulento) 


o — massa específica do fluido 


c,— calor específico do fluido a pressão 
constante 


u — velocidade 
S — secção do escoamento 
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Re — número de Reynolds 


P, = o * = número de Prandtl 
N,= me = número de Nusselt 


Pe == Re. Rr — número de Péchet 
= fluxo térmico por unidade de comprimento 
do tubo. 
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NOTICIÁRIO 


Colóquio Internacional 
Sistemas lógicos. Concepção e aplicações 


A Société Royale Belge des Electriciens organiza 
em Bruxelas nos Laboratórios de Electrónica 
Industrial e Automática da Universidade livre 
de Bruxelas de 15 a 20 de Setembro 1969, 
um colóquio internacional sobre SISTEMAS 
LÓGICOS, sua concepção por meios sistemáticos 
e suas aplicações científicas e industriais. 

Esta manifestação tem por fim criar contactos 
entre as pessoas que se interessam por um lado 
pelo estudo teórico, e por outro pela sua apli- 
cação científica e industrial. 

Informações : 


Dr. J. Florine 

Colloque International 

Université Libre de Bruxelles 50, avenue F. D. 
Roosevelt Bruxelles 5 — Belgique 


Da 


Conferências PRO.AQUA — «A água 
e o ar na indústria» 


Como nas precedentes realizações, em 1958, 
1960 e 1965 o 4.º salão internacional sobre ar e 
água, PRO-AQUA 69, que terá lugar de 29 de 
Maio a 4 de Junho de 1969, nos salões da Feira 
Suíça de amostras em Bale, será acompanhada 
de conferências internacionais. As conferências 
gravitarão à volta do tema «A água e o ar na 
indústria», e 3 temas distintos formarão 3 grupos: 

No tema A os especialistas Alemães, Franceses, 
Ingleses, Belgas e Suíços falarão do estado actual 
das realizações de instalações de purificação de 
águas usadas e da utilização de água em circuito 
fechado nas indústrias dos respectivos países. 

O tema B é consagrado ao saneamento do ar 
na indústria. 

O tema C dirá respeito aos problemas de dis- 
tribuição de água e da eliminação das águas 
usadas nalguns grupos industriais. 

De 2 a 5 de Junho de 1969 desenvolver-se-á, 
em relação com a PRO-AQUA 69, o 4.º congresso 
internacional sobre os tratamentos das imundií- 
cies (GIROM). 


V Colóquio de Directores e Técnicos de 
Fábricas de Cimento 


Durante os dias 5 a 8 de Maio próximos, 
realizar-se-ão, no Instituto Eduardo Torroja da 
Construção e do Cimento, os V Colóquios de 
Directores e Técnicos de Fábricas de Cimento, 

O tema geral deste V colóquios será a «Auto- 
mação de Fábricas de Cimento», e até ao momento 
verifica-se um grande interesse principalmente 
nos países de língua espanhola e portuguesa. 

Trabalha-se actualmente na organização dos 
colóquios que reunirão à volta de 200 especialistas 
de diversos países. 


Feira de Hannover de 1969 
26 de Abril a 4 de Maio 


Apresentam-se nesta feira (pavilhão 7) pro- 
dutos de alta qualidade em ferramentas, sendo 
a participação estrangeira de quase 1/3, perfa- 
zendo um total de 101 casas de 19 países. Esta 
cifra compõe-se de 27 expositores directos e 74 
casas representadas. Estas firmas provêm essen- 
cialmente dos Estados Unidos, França, Grã-Bre- 
tanha, Países Baixos, Suécia e Suiça. 

No pavilhão 5-A apresenta-se uma exposição 
sobre máquinas de ensaio e balanças. Expõem 
26 expositores directos e 8 firmas, pertencendo 
à Gran-Bretanha, Áustria, Suécia, Suíça e Esta- 
dos Unidos. 

Nos edifícios 11 e 11-B, 105 casas expõem 
aparelhos de Rádio e Televisão, assim como arti- 
gos de Fonotécnica e antenas. Além destes pavi- 
lhões, as indústrias eléctricas ocupam ainda os 
edifícios 1, 10, 11, 11-A, 12 e 13. Encontram-se 
aqui grupos elaboradores de produção e trans- 
formação de electricidade, distribuição e aplica- 
ção de energia eléctrica técnica de informação, 
técnica de recepção de rádio, acústica eléctrica, 
elementos de construção, técnicas de medição, 
controle, regulação e distribuição. 

Nas salas 16-A, 16-B e 16-C estão expostos 
artigos sobre Oleohidrosílica e Pneumática, bom- 
bas, condensadores, e armaduras. 


TECNICA N.º 385 
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C. D. U. 681.14:218.3 


CALCULADORES ANALÓGICOS 


APLICAÇÕES A SISTEMAS MECÂNICOS* 


RESUMO 


Este artigo trata da aplicação de calculadores analó- 
gicos a sistemas mecânicos, São analisados os seguintes 
uspectos : 

— Análogos eléctricos, 

— Tipos de análogus electricos. 

— Circuitos lineares de parâmetros concentrados, 

— «inálogo da relação entre forças e deslocamentos 
generalizados em elementos actuados em dois pontos, Aplica- 
ção a sistemas complexos, 

— Aspectos conceptuais e metodológicos do cúlculo 
analógico. 

— Anúálogos de alguns sistemas mecânicos simples em 
regime estático e dinâmico, 

— Análogus de uma viga em regime estático e dinâmico, 

— Análogo de uma barra sujeita a torção em regime 
dinâmico, 

— Regimes livres, Frequências próprias, Representação 
modal, Análise de vibrações. 

— Anúlogos de elementos não lineares, 

— Análogos de placas em regime estático. 

— Análogos aplicáveis a problemas de condução de calor. 


1. CÁLCULO ANALÓGICO 


O cálculo analógico corresponde na sua essên- 
cia a uma atitude mental, como qualquer activi- 
dade intelectual do homem. 

Dada a facilidade de aplicação, comodidade e 
baixo preço, tem sido objecto de um sem número 
de utilizações dispersas, nos mais variados domí- 
nios e com os mais diversos pontos de vista e 
meios. 

Consequentemente, nesse quase caos de formas 
concretas, sobressai com relativa nitidez a essência 
das ideias, no sentido de Platão. E, da generali- 
dade de assuntos tratados, sobressai a: grande 
uniformidade dos domínios da ciência e da téc- 
nica, quanto aos seus aspectos fundamentais. 

O cálculo analógico é a transposição para a 
ciência e técnica da filosofia das ideias de Platão, 
sendo surpreendente a profundidade da sua cor- 
respondência com essa filosofia, o que traduz, 


Por CARLOS DE MEDEIROS PORTELA 


Doutor em engenharia electrotécnica 


SYNOPSIS 


This paper deals with analogic computers and their 
application to mechanical systems. The following aspects 
are presented ; 

— Electrical analogs. 

— Types of electrical analoga. 

— Linear networks of lumped parameters, 

— Analoo of the relation between generalized forces and 
displacements in elements actuated in two points, Applica- 
tion to complex systems. 

— Conceptual and methodologic aspects of analogie 
computation, 

— Analogs of some simple mechanical systems in static 
and dynamic regime. 

— Analogs of a beam in static and dynamic regime. 

— Analogs of a beam submited to torsion in dynamic 
regime, 

— Free regimes. Natural frequencies. Modal repre- 
sentation, Vibration analysis. 

— Analogs of non linear elements. 

— Analogs of plates in static regime, 

— Analogs applicable to heat conduction problems. 


talvez, a causa da fecundidade do cálculo analó- 
gico, principalmente como atitude mental, para 
além das formas particulares de que se revista a 
sua aplicação. 

E, assim, não é de admirar que se chegue a 
beneficiar de grande parte das vantagens do cál- 
culo analógico utilizando apenas as suas atitudes 
e métodos topológicos, matemáticos e lógicos, 
mesmo sem dispor de análogos físicos concretos. 
Por exemplo, o cálculo analógico pode ser um 
simples processo de formular um problema para 
a sua resolução num calculador digital. 

Quando, para determinar a altura das pirâmides 
do Egipto, um grego utilizou uma estaca e a 
medida de comprimento da sua sombra, estava a 
aplicar uma forma já relativamente evoluída do 
cálculo analógico, mais abstracta que a simples 
escala geométrica, pois, além dessa escala, utilizou 
uma estilização dos parâmetros da pirâmide que 
condicionam a sua sombra. É curioso que nessa 


(*) Conferência proferida no Instituto Superior Técnico em 10 de Fevereiro de 1969. 
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aplicação concreta tenha surgido a sombra física, 
tal como na alegoria platónica. 

No seu sentido mais comum, pode dizer-se 
que se utiliza o cálculo analógico quando se esta- 
belece uma correspondência entre os parâmetros 
de um fenómeno a estudar e os parâmetros de 
um sistema físico concretamente representado. 

Deve notar-se que, abstraindo um pouco, se 
pode prescindir da condição de representar con- 
cretamente o sistema físico análogo, sem perder 
o que há de essencial e mais útil no cálculo 
analógico. 

O fenómeno físico escolhido para análogo pode 
ser qualquer. A régua de cálculo é um exemplo 
simples de calculador analógico. Estabelece-se 
uma correspondência entre dois números e os 
seus logaritmos, e efectua-se mecânicamente a 
soma dos logaritmos, «análogo» da multiplicação 
dos dois números. 

Verificou-se um desenvolvimento maior e mais 
sistemático de análogos eléctricos, basicamente 
por razões tecnológicas, de custo, e de facilidade 
de variação de parâmetros e de estabelecimento 
de ligações topológicas. 

No entanto, a escolha de grandezas electro- 
magnéticas para o análogo nada tem de essencial. 
Essa escolha, se por um lado foi e é vantajosa, 
teve, no entanto, um inconveniente, ligado ao 
carácter degenerado de algumas relações tenso- 
riais dos circuitos eléctricos lineares estáticos 
passivos. Essa degenerescência tem diversas van- 
tagens quanto à simplicidade das relações, mas, 
em certa medida, constituiu um obstáculo ao 
reconhecimento do carácter intrinsecamente ten- 
sorial de algumas relações, e constitui uma difi- 
culdade na representação física com circuitos 
eléctricos concretos e simples de relações em que 
se não verifique essa degenerescência. 


2. TIPOS DE ANÁLOGOS ELÉCTRICOS 


Para muitas aplicações concretas, a utilização 
de um análogo particular muito simples é a 
forma mais adequada de resolver problemas por 
vezes complexos. 

Este facto, se por um lado é útil, tem sido o 
maior obstáculo ao desenvolvimento sistemático 
e generalizado do cálculo analógico. A situação 
quase oposta dos calculadores digitais constituiu 
a causa do relativo atraso com que se desenvol- 
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veram, mas constitui também a causa da sua 
rápida generalização, uma vez atingido o «limiar» 
de justificação que permita a sua existência. 
Há grande diversidade de análogos eléctricos, 
sendo um tanto artificial a sua sistematização. 
Podem considerar-se, basicamente, os seguintes 
grandes grupos : 
— Circuitos lineares de parâmetros concentra- 
dos. 
— Circuitos não lineares de parâmetros con- 
centrados. 
— Circuitos de parâmetros distribuídos. 
— Circuitos com relações algébricas num 
domínio complexo. 
— Circuitos com elementos normalmente desig- 
nados por amplificadores operacionais. 
Em palestra anterior foram apresentadas as suas 
principais características (!). Por esse motivo, e 
para não alongar a exposição não o faço aqui. 


3. CIRCUITOS LINEARES DE PARÂMETROS 
CONCENTRADOS 


Vou indicar em primeiro lugar algumas ana- 
logias eléctricas aplicáveis ao caso de estruturas 
mecânicas, utilizando análogos constituídos por 
elementos lineares passivos de parâmetros con- 
centrados. 

Na fig. 3-1 indicam-se algumas relações para 
elementos simples do tipo anterior. 


4. ANÁLOGO DA RELAÇÃO ENTRE FOR- 
CAS E DESLOCAMENTOS GENERALIZA- 
DOS EM ELEMENTOS ACTUADOS EM 
DOIS PONTOS. APLICAÇÃO A SISTEMAS 
COMPLEXOS 


Consideremos em primeiro lugar, para realçar 
alguns aspectos essenciais do cálculo analógico, 
o caso de tubagens hiperestáticas, por exemplo 
como representado esquemáticamente na fig. 4-1, 
constituídas por três ramos ab”, a” b”,a” b”, 
designados genéricamente por ab. 


Cada ramo é caracterizado, em determinadas 
condições, por uma relação da forma 


1) 4X — Xp |=1Th.FI 


(') «Introdução ao cálculo analógico”. Técnica 376, 1968. 
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AROS E Li 


E qr é SU 
Fr 
Fig. 4-1 — Tubagens hiperestáticas a' b',a” b" ea” b” 
ligando o corpo B, e a união A. As três tubagens são 
actuadas nos pontos a', a”, a”' por forças generalizadas 


(forças e momentos) definidas pelas matrizes | F' ||, 
EC, EPI. 


sendo 
IX matriz de coordenadas generalizadas 
do extremo do ramo, em geral cons- 
tituída por seis elementos, sendo três 
as coordenadas cartesianas e três as 
rotações segundo o mesmo sistema de 
coordenadas; 
matriz de forças generalizadas actuan- 
do num extremo do ramo, em geral 
constituída por seis elementos, sendo 
três as componentes cartesianas da 
força e três as componentes do mo- 
mento, 


IF |) 


Ter-se-á, ainda, 


4) IF|l=|S1. 


[4XI, 
sendo || S || a matriz inversa da matriz | T ||. 


Procuremos em primeiro lugar o análogo de 
cada ramo, considerado isoladamente. Pode ter-se, 
por exemplo, a representação indicada na fig. 4-2, 
fazendo corresponder as forças e deslocamentos 
a correntes, em regime alternado sinusoidal, 
estando em fase as grandezas do mesmo sinal e 
em oposição as grandezas de sinal oposto. 

Nota-se que, estando todas as correntes alter- 
nadas em fase ou oposição, 2s correntes podem 
ser tratadas e medidas como grandezas algébricas 


constantes, abstraindo do seu carácter sinusoidal 
no tempo (!). 

Por comodidade de notação, as correntes aná- 
logas de AX, e de F são representadas por 4X, 
e E,, admitindo, naturalmente, factores de «escala» 
adequados. 

Este circuito constitui o análogo de um ramo 
e impõe entre as doze variáveis das matrizes 
|4XI|| e IF] a relação 2). À parte essa rela- 
ção, as doze grandezas são livres. Salienta-se que, 
no análogo, tal como no sistema, há seis relações 
entre doze grandezas e seis graus de liberdade. 

Vejamos, por exemplo, como, a partir deste 
análogo, se podem determinar imediatamente os 
36 elementos da matriz || S ||, isto é, inverter a 
matriz || T||. Impondo a corrente 


6) AX,=1, 


7) AX=...=AX=0, 


e tornando livres as 6 «correntes» F,,. .. E,, 
curto-circuitando os seis pares de terminais cor- 
respondentes, os valores de F,,. .., F; nestas 
condições são iguais aos elementos de uma 
coluna da matriz inversa || S. || Fazendo, sucessiva 


e independentemente, igual a 1, AX,,..., AX pá 


com as restantes variáveis À X, nulas, obtém-se, 
por medida directa de corrente, os 36 elementos 
da matriz inversa. 

Se, por exemplo, forem impostas três das 
variáveis F. e três das variáveis X,, basta injectar 
as correntes respectivas nos terminais correspon- 
dentes, e curto-circuitar os pares de terminais 
relativos às outras seis variáveis. 

O anélogo representa as relações entre as doze 
variáveis impostas pelo ramo. Essas doze variá- 
veis são doze correntes em pares de terminais, 


(1) Supondo, por exemplo, todas as grandezas, y;, da 
forma 


5) y=Y, cos (vt + 4), 


sendo t o tempo e w e v os mesmos para todas as gran- 
dezas y;, o parâmetro Y, define completamente y, , poden- 
do constituir o análogo de uma grandeza de carácter algé- 
brico, independente de t. Designa-se, nesta exposição Y,; 
por «amplitude algébrica» da grandeza sinusoidal y,, num 
circuito em corrente alternada em que as tensões e correntes 
sejam todas da forma 5), estando portanto em fase as do 
mesmo sinal e em oposição as de sinal contrário. 
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| Transformador 1 


Fig. 4-2 — Análogo de um ramo genérico ab das tubagens representados na 
fig. 4-1. O análogo é constituído por seis transformadores (ideais), 1 a 6, 
cada um de sete enrolamentos. Os valores 1, T,,, Tj9..., indicados junto 
de cada enrolamento, são valores proporcionais aos parâmetros n, desses 
enrolamentos. As correntes F,,..., E, são sinusoidais, de pulsação w, em 
fase ou oposição (anulando-se todas simultâneamente) de amplitudes (em 
sentido algébrico) proporcionais às componentes F,,,..., F; da matriz || F | 
das forças generalizadas que actuam o ramo no ponto a. As correntes 
AX,,...s AX; são correntes sinusoidais de pulsação «w (anulando-se 
simultâneamente com as correntes F,) de amplitudes (em sentido algébrico) 
proporcionais às componentes AX, ..., AX, da matriz || AX || das deforma- 
ções generalizadas (diferença entre os deslocamentos generalizados dos 
pontos a e b) do ramo ab. 


que, à parte as relações impostas pelo ramo, são 
livres. Portanto, a constituição interna do aná- 
logo do ramo pode ser ignorada ao compor as 
condições externas de ligação do ramo ao exterior. 

Para essa composição há, agora, que simular 
as ligações exteriores, de maneira «análoga» às 
ligações entre os ramos e as forças generalizadas 
exteriores que os actuam. 

Os vários conjuntos de doze pares de terminais 
do circuito correspondentes a cada ramo devem 
ser ligados por forma a impor as relações entre 
os deslocamentos generalizados e forças genera- 
lizadas impostas «externamente» pela forma como 
os ramos estão ligados. 

Suponhamos, por exemplo, que é imposta 
externamente a componente, segundo uma das 
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direcções (relativa ao índice 1 dos elementos das 
matrizes ||Fl||) da resultante das forças que 
actuam as extremidades a das barras, isto é, que 
é imposto externamente o valor 


8) E =F + F, + Me 


A «soma» das componentes F, das forças nas 
três barras, tem por análogo a soma das correntes 


hd t tel 
análogas de E, EF. F,- Para obter essa soma, 


basta ligar em paralelo os três pares de terminais 
correspondentes respectivamente a essas corren- 
tes e impor, por meio de um gerador de corrente, 


que essa soma tenha o valor «análogo» de E p 
conforme representado na fig. 4-3. 
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Fig. 4-3 — Análogo do conjunto de tubagens hiperestáticas constituídas pelos ramos a'b', a” b”, a” b” da fig. 4-1. 
Representação das condições de ligação entre os três ramos. O análogo de cada ramo é considerado, para efeito 
destas condições de ligação, apenas por meio dos doze pares de terminais correspondentes às seis componentes 
da matriz de deformações generalizadas || AX || e às seis componentes da matriz das forças generalizadas || F ||. 
Os geradores de corrente relativos às condições impostas não nulas (neste exemplo Fr, A Xo, P", e P”,) são 
geradores de corrente alternada da mesma pulsação «w e tais que as várias correntes sejam simultâneamente 
nulas, estando em fase ou oposição (consoante o sentido das grandezas que representam). Por medida das 
correntes respectivas, obtém-se directamente os elementos das matrizes || AX || e I||Fi| relativas aos vários 


ramos (considerando, naturalmente, os factores de escala adoptados). 
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Se for imposta a condição de as três compo- 
nentes da deformação das barras segundo a 
direcção 1 serem iguais, 


9) AX, = 4X =AX, , 


basta impor que as três correntes «análogas» a 
Â X, ; X,» A, X, sejam iguais, para o que basta 
ligar em série os três pares de terminais corres- 
pondentes, como representado na fig. 4-3. 

Na fig. 4-3 representam-se ainda, a título 
exemplificativo, as ligações relativas às seguintes 
hipóteses : 

a) Componentes segundo uma direcção (repre- 
sentada pelo índice 2) das forças F, F”, PF”, 
livres. 

b) Deformações segundo uma direcção (repre- 
sentada pelo índice 2) iguais entre si e iguais a 


H i 
um valor imposto externamente, À X ge 


c) Componente segundo uma direcção (repre- 
sentada pelo índice 3) da resultante das forças 
F, F”, PF” iguais a zero. 

d) Diferença imposta (por meio de um pré- 
-esforço, por exemplo) entre as deformações dos 
três ramos, segundo a direcção correspondente 
ao índice 3. 

e) Momento resultante dos momentos actuando 
nos três ramos (representados pelos elementos de 
índices 4; 5 e 6 das matrizes || F||) nulo. 

f) Rotações dos três ramos (representadas 
pelos elementos de índices 4; 5 e 6 das matrizes 
| 3X ||) iguais. 


5. ASPECTOS METODOLÓGICOS DO 
CÁLCULO ANALÓGICO 


Há os seguintes aspectos do estabelecimento 
do análogo que convém realçar : 


1) Decomposição do sistema físico em partes, 
admitindo cada uma. ligada às restantes 
por variáveis de ligação definidas, neste 
caso as deformações dos ramos e as forças 
num nó de ligação, podendo essas variáveis 


ter valores quaisquer. 
2) Consideração de cada uma destas partes 


como um sistema físico independente, 
abstraindo das relações entre as variáveis 
de ligação impostas pelo resto do sistema. 

3) Estabelecimento do análogo de cada uma 
das partes isoladamente. 
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4) Estabelecimento do análogo das relações 
entre as variáveis de ligação correspon- 
dentes às ligações entre as várias partes, 
de maneira quase meramente topológica, 
abstraindo dos condicionamentos impostos 
por cada uma das partes do sistema. 


Em casos simples, estes vários aspectos apare- 
cem misturados, não havendo interesse especial 
em separá-os. 

Em sistemas complexos, no entanto, este método 
de «corte» do sistema em subsistemas elemen- 
tares, e consideração separada do comportamento 
de cada elemento, por um lado, e das relações 
entre variáveis de ligação, por outro lado, é extre- 
mamente importante. O problema das relações 
entre variáveis de ligação tem em grande parte 
um carácter topológico que, se devidamente uti- 
lizado, se adapta muito bem ao análogo eléctrico. 
A ligação de dois terminais, por exemplo, é uma 
forma extremamente simples de impor uma con- 
dição topológica. 

No que respeita básicamente aos circuitos eléc- 
tricos, o problema geral deste tipo de métodos, 
sob os pontos de vista filosófico, físico, topoló- 
gico, geométrico e matemático, tem sido objecto 
de análise nos últimos tempos. 

Especialmente sob o ponto de vista das aplica- 
ções à electricidade, há uma exposição relativa - 
mente sistemática de Kron, que usa para c 
método a palavra «diakoptics», que pretende 
significar «sistema de cortar». 

Não vem ao caso discutir aqui as teorias de 
Kron; há aliás algumas das suas afirmações de 
que discordo. De qualquer forma, no entanto é 
extremamente importante a unificação, ou pelo 
menos conjugação, para o domínio da física e da 
técnica, deste tipo de atitude, ou método cientií- 
fico, que, aliás, está de acordo com uma das regras 
do método de Descartes: dividir a dificuldade 
em tantas partes quantas as necessárias para 
separadamente melhor as resolver. 

Esta unificação ou conjugação constitui um 
tema que está em grande parte inexplorado, talvez 
mais em engenharia mecânica que electrotécnica. 


6. ANÁLOGOS DE ALGUNS SISTEMAS 
MECÂNICOS SIMPLES EM REGIME ESTÁ- 
TICO E DINÂMICO 


Apresentado o exemplo analisado em 4., com 
relativa generalidade, passo a indicar alguns 
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exemplos muito simples de análogos eléctricos 
de sistemas mecânicos. 


a) Mola linear unidireccional 


Consideremos uma mola linear unidireccional, 
caracterizada, mum espaço a uma dimensão, x, 
por dois pontos 1 e 2 de coordenadas X, e X,, 
por forças F, segundo a direcção x, actuando em 
sentidos opostos nos pontos 1 e 2, e pela relação 


1) F=K(X—X), 


sendo K, constante, uma característica da mola. 

Na fig. 6-1 representa-se um análogo eléctrico, 
adequado para regime estático da mola, sendo 
1 e 2 os terminais de uma resistência R, e U, 
e U, as tensões dos terminais 1 e 2 em relação 
a um terminal de referência de tensões, O. Em 
regime estacionário ou em regime alternado 
simusoidal e com todas as correntes e tensões em 
fase ou oposição, tem-se a relação 


2) 1=)(W,-U), 
R a 


sendo, no caso de grandezas sinusoidais, Ie U 
as «amplitudes algébricas» das tensões e correntes. 
Ter-se-á, portanto, a seguinte correspondência 
entre as grandezas mecânicas e eléctricas, para 
o análogo anterior (utilizável, em sistemas 
complexos, em regime mecânico estático): 


Grandeza mecânica Grandeza eléctrica 


— 
DDD >> 


deslocamento X tensão U 
força F corrente 1 
K a 
R 
a . 1 
energia elástica potência vezes — 


Em regime dinâmico, para a sua integração 
num sistema mecânico complexo, o análogo 
anterior não é, em geral, adequado, sendo pre- 
feríveis outras analogias, uma das quais represen- 
tada na fig. 6-2. 

Sejam v, e v, as velocidades dos pontos 1 e 2 
da mola (segundo a direcção x). Tem-se a relação 


3) E=K [ (uv) dt. 


No circuito eléctrico da fig. 6-2 tem-se a 
relação 


4) À == Tu) dt. 


Se, na mola, houver uma dissipação de tipo 
viscoso, caracterizada por uma potência de dissi- 
pação da forma 


Mola linear unidireccional 


Simbologia 
eee 
Relação entre as 
grandezas ER pe] 
X 
Correspondência F 
entre as grandezas K 


Análogo eléctrico 


R 


+ dA + 
a ; “4, 
Ss E 
É) 


Fig. 6-1 — Análogo eléctrico de uma mola linear unidireccional em regime estático. As grandezas eléctricas 

(tensões U,, U, e corrente 1) podem ser estacionárias, amplitudes (em sentido algébrico) de grandezas alternadas 

sinusoidais, todas da mesma frequência e em fase (ou oposição se de sinais contrários), ou valores instantâneos 
num mesmo instante. 
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Mola linear unidireccional Análogo eléctrico 
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Fig. 6-2 — Análogo eléctrico de uma mola linear unidireccional em regime dinâmico. A resistência R existirá 
apenas no caso de se considerar um efeito dissipativo na mola, de tipo viscoso. No caso contrário, 


são nulas 


5) Pa = B(m—v)”, 
e por uma força de dissipação da forma 


6) f=B(v,—v,), 


essa dissipação será correctamente representada 
no análogo por uma resistência R, tendo-se uma 
potência de dissipação 


1 
7) p= " (u, — u)) A 


e uma corrente (análogo da força) da forma 


4 


1 
8) E ia 


Portanto, ter-se-á a seguinte correspondência 
entre as grandezas mecânicas e eléctricas, no aná- 
logo dinâmico anterior : 
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f, i', B. 


Grandeza mecânica | Grandeza eléctrica 


velocidade v tensão u 
força E corrente i 
1 
K guenísio 
L 
integral da tensão 
deslocamento X . 
[udt 
tempo t | tempo t 
potência potência 
energia | energia 
B 2 
DR 


b) Barra rígida actuada por forças paralelas 


Suponhamos uma barra rígida actuada por 
forças paralelas, conforme representado nas figs. 
6-3 e 6-4, sendo X,,...,X, as coordenadas x 
de n pontos da barra segundo a direcção do deslo- 
camento (aproximademente ortogonal à barra), 
Vj-«., V, às velocidades correspondentes, M,,..., 
M, as massas de troços consecutivos da barra na 
vizinhança dos n pontos e supostas concentradas 
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neles, e F,..., F, forças paralelas à direcção 
x actuando nos n pontos. 

Suponhamos em primeiro lugar a barra em 
equilíbrio estático. Ter-se-á, com a notação da 
fig. 6-3, 


A, — X, ] 

9) E Ls itis=taia=s) 
PO, | 
i—1 


10) LEFMAE + + LI E+ 
k = 1 
n—1 


a ad e e ( à Lo E O 


k = | 


11) E+E+H...tE+...+E=o0. 


Consideramos um transformador ideal de n-1 
enrolamentos, ligados em série, conforme se 
indica na fig. 6-3, e a notação das tensões e cor- 
rentes indicada na mesma figura. Suponhamos 
todas as correntes e tensões alternadas sinusoi- 
dais, de pulsação », anulando-se simultâneamente 
todas as grandezas, que estarão em fase se do 
mesmo sinal e em oposição se de sinal contrário. 

Sejam U, e I, as amplitudes (em sentido algé- 
brico) das tensões e correntes indicadas. Ter-se-á 


12) - —=— 


| —1 
dd mito Ends dd E ndi+ 
k = 1 


++ n)h=0 


k=1 
14) Elo Fi. ccti=0: 
Pode, assim, estabelecer-se a seguinte corres- 
pondência entre as grandezas do sistema mecà- 


nico, em equilíbrio estático, e as grandezas do 
análogo eléctrico (em regime alternado sinusoidal): 


Grandeza eléctrica | 


tensão U 


| Grandeza mecânica 


deslocamento X 
| força F 
| 


corrente | 


relação de braços relação de transfor- 


I, 
| 
| n, 
mação —— 
ni 
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Suponhamos agora a barra em regime dinâmico, 
e seja v, a velocidade de um ponto genérico i 
(segundo a direcção x), 


15) v, a 
dt 


Ter-se-á, com a notação da fig. 6-4, 


X, — X, 1. 
| EL es (,j=1;2;...;n—1) 
X 41X , 
dv, 
17) F=M,-—— (i=1:2: sn) 


18) IME) +HML+HLME-E)+...+ 


, 
ae 


FOZIJE-EI E... + 


k = | 


n—1l 


+(ZIN+H(E-E)=0. 


k = | 


19) F-ER+HE-E+...+tE—F=o0. 


Consideremos o análogo eléctrico representado 
na fig. 6-4, constituído por um transformador de 
n-1 enrolamentos e n capacidades C,,..., € 
Tem-se, com a notação da mesma figura, 


nº 


n 
9») E = (i, j=1;...;n-1) 
n 


21) ;=C —— (,j=1;...;n=—1) 


22) ni) ++) (i++ 
i— 1 n— 

+ nd) +... HS ny (,-i)=0 
k=1 k=1 


2) Et sti,—i=0 


Portanto, pode estabelecer-se a seguinte corres- 
pondência entre as grandezas do sistema mecà- 
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nico em regime dinâmico e as grandezas do 
análogo eléctrico : 


Grandeza mecânica | Grandeza eléctrica 


velocidade v tensão u 
deslocamento x f udt 
tempo t tempo t 
força F corrente i 
massa M capacidade € 
força F' corrente i' 
aceleração add sa 

ma 7 dt 


relação de braços 


L, 


relação de transforma- 


a ção -e 

I; n | 
potência potência 
energia | energia 


c) Combinação dos exemplos a e b 


Quer no caso de regime estático quer no de 
regime dinâmico, a correspondência entre as gran- 


dezas mecânicas e eléctricas nos exemplos a e b 
é a mesma. Consequentemente, podem obter-se 
análogos de sistemas constituídos por combinação 
de barras rígidas e molas, compondo os respecti- 
vos análogos e impondo as relações adequadas 
entre as variáveis de ligação respectivas. 

Seja, por exemplo o caso de uma barra actuada 
por três molas e por uma força exterior F no 
ponto 1, conforme representado esquemáticamente 
nas figs. 6-5 e 6-6. Têm-se, respectivamente em 
regime estático e dinâmico, os análogos indicados 
nessas duas figuras. No caso estático, a força 
exterior F é simulada por um gerador de corrente 
alternada siunsoidal, de amplitude complexa F. 
No caso dinâmico, F é simulada por um gerador 
de corrente variando no tempo de maneira aná- 
loga a F (à parte o factor de escala dos tempos 
do sistema mecânico e do análogo eléctrico e o 
factor de escala de conversão de forças no sistema 
mecânico em correntes no análogo). 


7. ANÁLOGOS DE UMA VIGA EM REGIME 
ESTÁTICO E DINÂMICO 


Consideremos uma viga actuada por forças 
unidireccionais paralelas entre si e perpendi- 


Sistema mecânico 
(em regime estático) 


Simbologia 


Correspondência 
entre as 
grandezas 


Análogo eléctrico 


Fig. 6-5 — Análogo de um sistema constituído por uma barra rígida e três molas, actuando por uma força, em 
regime estático. 
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Sistema mecânico 


A Análogo eléctrico 
| (em regime dinâmico) 


Simbologia na 
a | n 
d, Do 
| 1 
. | K, , Es: 
Correspondên- i 
cia entre as M, C, 
grandezas 4. u, 
] 
X; | fu dt 
t t 
F i 
Fig. 6-6 — Análogo de um sistema constituído por uma barra rígida e três molas, actuado 
por uma força, em regime dinâmico. 
culares à viga. Para concretizar, suponhamos o Consideremos uma estilização da viga, assimi- 
caso representado na fig. 7-1a. lando-a a um conjunto de peças rígidas em 


" 
a) ema 


K = EL M-mãÃx v=[5— 


(26 
Sc = [3 


K.const, 


b) 


Fig. 7-1 a) Viga actuada por forças unidireccionais paralelas entre si e perpen- 
diculares à viga. 
b) Estilização da viga, assimilando-a a um conjunto de peças rígidas 
em forma de T, articuladas entre elas, ligadas por molas, supondo 
concentradas nas articulações as massas de troços da viga. 
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forma de T, articuladas entre elas, e ligadas por 
molas, supondo concentradas nas articulações 
as massas M de troços da viga correspondentes 
a um invervalo 4x, conforme representado na 
fig. 7-1b. Esta estilização corresponde, como se 
verifica facilmente, a uma aproximação por dife- 
renças finitas das equações diferenciais da viga, 
para pequenas deformações. 

Na fig. 7-2 e 7-3 representam-se os análogos 
do comportamento de uma das peças rígidas em 


Fig. 7-2 — Análogo do comportamento de uma das peças 
rígidas em T da estilização da viga indicada na fig. 7-1 b, 
em regime estático. Os pontos a e b são as articulações 
do T aos dois T consecutivos e c é a ligação às molas 
que ligam o T aos dois T consecutivos. S, e Sy repre- 
sentam as forças transversais nas secções ae b da viga 
referidas ao troço a be F,. 1 (sendo 1/a unidade de 
comprimento) representa a momento flector na secção 
c (intermédia entre as secções a e b). Os deslocamentos 
dos pontos a e b segundo Y são representados por 
Ta e Yge 

O análogo é constituído por um transformador ideal 
de dois enrolamentos, de relação de transformação 1: 
Ax, em regime alternado sinusoidal, estando todas as 
tensões e correntes em fase, ou oposição, consoante 
os seus sinais relativos. As grandezas S,, Sp, Fo têm 
como análogos as amplitudes algébricas das correntes 
Sa Sp: Ep € as grandezas Y,, Yp, 9 as amplitudes 
algébricas das tensões Ya, Y,, 9, (em relação a um nó 
de referência, O, no que respeitaa Y, e Y,). 

O análogo traduz as relações: 


Ya— Yp 
Pa SRS A 
AX ç 
a -— Sp 
Fo Sa + Sp AX 
2 


T, em regime estático e dinâmico, constituído 
um transformador de dois enrolamentos. Na 
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Fig. 7-3 — Análogo do comportamento de uma das peças 
rígidas em T da estilização da viga indicada na fig. 7-1 b, 
em regime dinâmico. Os pontos a e b são as articula- 
ções do T aos dois T consecutivos e c é a ligação às 
molas que ligam o T aos dois T consecutivos; S, e Sp 
representam as forças transversais nas secções a e b 
da viga referidas no troço ab e F,. 1 (sendo 1a unidade 
de comprimento) representa o momento flector na 
secção c (intermédia entre as secções a e b). Às velo- 
cidades dos pontos a e b segundo Y são representadas 
por Y,,€ Ysp: 

O análogo é constituído por um transformador ideal 
de dois enrolamentos de relação de transformação 
r:4x, em regime não estacionário. As grandezas 
Sa Sp Fe. têm como análogo os valores instantâneos 
das correntes S,, Sp, Fo € as grandezas V,, Vpy Op OS 
valores instantâneos das tensões v,, Vp, vç (em relação 
a um nó de referência, O, no que respeita a v, e v,). 
O análogo traduz as relações 


fig. 7-4 representa-se o análogo do comporta- 
mento de uma mola da estilização da viga, em 
regime estático, constituído por uma resistência, 
e na fig. 7-5 o análogo em regime dinâmico, 
constituído por uma indução. Na fig. 7-6 repre- 
senta-se o análogo do comportamento da massa 
de um troço de viga, suposta concentrada na 
articulação de duas peças rígidas em T, em 
regime dinâmico. Para facilidade de notação, as 
tensões, correntes, e relações de transformação 
são representadas com a mesma notação das 
grandezas mecânicas de que são análogas. 

Na fig. 7-7 representa-se o análogo da viga, 
em regime estático. As forças F, e F, são simu- 
ladas por dois geradores de corrente alternada 
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Sistema mecanico 


Analogo electrico 


Fig. 7-4 — Análogo do comportamento de uma mola da 
estilização da viga indicada na fig. 7-1b, em regime 
estático. o, e Oy representam os valores de & nas 
secções c e d da viga e F.l (sendo 1 a unidade de 
comprimento) o momento fletor no ponto a, intermé- 
dio entre os pontos c e d, O análogo é constituído por 
uma resistência em regime alternado sinusoidal, estando 
todas as tensões e correntes em fase, ou oposição, 
consoante os seus sinais relativos. As grandezas O, e 
Og têm como análogos as amplitudes algébricas das 
tensões q e 94, em relação a um nó de referência, 6, 
é a grandeza F (ou F.1) a amplitude algébrica da cor- 
rente F (ou F,1). O análogo traduz a relação 


Fi=kK. (0,—0,). 


Sistema mecânico Anâlogo electrico 


Fig. 7-5 — Análogo do comportamento de uma mola da 
estilização da viga indicada na fig. 7-1b, em regime 


dinâmico w, e oq representam os valores de w, € my 
nas secções ce d da viga e F.I (sendo 1 a unidade de 
comprimento) o momento flector no ponto a, intermé- 
dio entre os pontos c e d. O análogo é constituído por 
uma indução. As grandezas q, e vq têm como análogos 
os valores instantâneos das tensões w, € 4, em reladão 
a um nó de referência, O, e a grandeza Il (ou F.l) o 
valor instantâneo da corrente F (ou F.1). O análogo 
traduz a relação 


p 
PL1I=K | (q — te ) dt. 


e 


sinusoidal (anulando-se simultâneamente), de 
amplitudes algébricas F, e F,. A condição de 
encastramento no apoio é simulada impondo 
o anulamento das tensões análogas de Y e & no 
ponto de encastramento. 
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Sistema Andalogo 


mecânico electrico 


Fig. 7-6 — Análogo do comportamento da massa de 
um troço da viga segundo a estilização indicada na 
fig. 7-1b, em regime dinâmico. A grandeza v representa 
a velocidade da massa M, suposta concentrada na 
secção a, segundo u direcção Y, e F a força de acele- 
ração correspondente. No análogo v é representada 
pela tensão v e F pela corrente F. 

O análogo traduz a relação 


e 


v 


F=M a 


Em regime estático F é nula, não sendo necessário 
incluir as massas M no análogo, 


Tem-se a seguinte correspondência entre as 
grandezas do sistema mecânico e do análogo 
eléctrico : 


Grandeza mecânica Grandeza eléctrica 


deslocamento Y tensão Y 
deslocamento angular & tensão & 
co É a 

AX R 
força F corrente I 
força transversal 5 corrente 5 


Na fig. 7-8 representa-se o análogo da viga, 
em regime dinâmico. As forças F, e F, são simu- 
ladas por geradores de corrente tais que a cor- 
rente varie em função do tempo de maneira 
análoga a F, e F,. 

Tem-se a seguinte correspondência entre as 
grandezas do sistema mecânico e do análogo 
eléctrico (Ver Quadro na página seguinte). 

No caso de haver dissipação de tipo «viscoso» 
dependente da velocidade v, essa dissipação pode 
ser representada no análogo por resistências em 
paralelo com as capacidades. Se houver dissi- 
pação de tipo viscoso dependente da velocidade 
angular w, essa dissipação pode ser representada 
no análogo por resistências em paralelo com as 
induções. 
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Fig. 7-7 — Análogo da viga representada na fig. 7-1, em regime estático. 
Este análogo resulta da composição dos análogos dos vários elementos 
da estilização do comportamento da viga indicados na fig.7-1b. As forças 
exteriores F, e F, são simuladas por geradores de corrente alternada 
sinusoidal (de igual frequência e anulando-se simultâneamente) gerando 
correntes de amplitudes algébricas F, e I,. A condição de encastra- 


mento é simulada impondo o anulamento de 9 e Y nos pontos corres- 
pondentes à secção de encastramento. 


Fig. 7-8 — Análogo da viga representada na fig. 7-1, em regime dinâmico. 

Este análogo resulta da composição dos análogos dos vários elementos da 

estilização do comportamento da viga indicados na fig. 7-1b. As forças 

exteriores F, e F, são simuladas por geradores de corrente tais que as 

correntes respectivas variem no tempo de maneira análoga a F,eF,. 

A condição de encastramento é simulada impondo o anulamento de v e q 
nos pontos correspondentes à secção de encastramento, 


8. ANÁLOGO DE UMA BARRA SUBMETIDA 


Grandeza mecêntes Grandeza eléctrica A TORC ÃO 
velocidade v = hi. tensão v 
dt foda Consideremos uma barra submetida a torção, 
sus v “ Als 
deslocamento Y = / vdt em regime dinâmico, conforme representado 
velocidade angular o = — tensão « esquematicamente na fig. 8-1. 
t A equação de deformação é 
deslocamento angular & J o dt poa ç 
| 
3x à [., de 
A DM(xgy) +—1!|K6)— |=I dia 
relação ca de braços do relação " transfor- ) (x) ds | (x) ) À (x) É , 
T rígido asia 
massa M =m 3x capácidade € 
nar ipi sendo aa deformação angular, t o tempo, e K 
força transversal 5 cousa um parâmetro de rigidez de torção, Io momento 
tempo tempo de inércia e M o momento de torção exterior, 
potência potência por unidade de comprimento (no que respeita a 
energia energia K Ie M) 
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' had ] e : 
1 1 

a — .— el 

n-1 nm a og 


Fig. 8-1 — Representação de uma barra submetida a torção 


A esta equação diferencial corresponde, sob a 


forma de diferenças finitas, a equação 


Kyo + 
2) M, dx + - (2541 -2) + 
À x 
Ko ' 
Y o n 
Cê AX (a, | %n) = 1, AX E Í 
sendo 
1;2;...;n-1;n;n+1;... 


indices identificadores de uma sucessão de secções 


1 “ A . 1 + - “ . 
da barra à distância 4x e n + 3 índice identi- 


ficativo de um ponto intermédio entre as secções 
nen+ti. 

Sob outra forma, naturalmente equivalente à 
anterior, pode considerar-se uma estilização da 
barra, representada na fig. 8-2, supondo-a cons- 
tituída por discos à distância Ax nos quais se 
suponha concentrada a inércia de um troço Ax 
da barra (1, 4 x), sendo dois discos consecutivos, 
relativos às secções n e n-+1, ligados entre si 
por uma mola de torção, caracterizada por um 

K 1 


n+ 


factor Rs de proporcionalidade entre a 


diferença de deslocamentos angulares nos extre- 

mos da mola e o binário de torção da mola. 
Consideremos o análogo representado na 

fig. 8-3, constituído por induções em série e 


Fig. 8-3 — Análogo de uma barra submetida a torção, em 
regime dinâmico. As correntes nv, obtidas por meio de 
geradores de corrente, são os análogos de 


fuçsa dt ; 


sendo M, o momento de torção exterior que actua um 
troço sx da barra na vizinhança do ponto n. À tensão u, 
é o análogo do deslocamento angular, «, na 
secção n da barra. 


capacidades em derivação, com a notação das 
grandezas aí indicada. 


Tem-se 


3) Pad li ! + au nm ql , 


donde 
do di, dt" di” 
4 PESE; a. 
j dt + dt + dt dt 
Ora, 
5) HE qo od ( Un) 
dt ac L E Upti— n 


K. o, 
Ca ? o E To Maca? 
Q 4 dimento e 
Ax , Lx Fiel 
A p . 
M. (0 K Mo! X Ma. Áx 


Fig. 8-2 — Estilização de uma barra submetida a torção. 
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6) dt - Lo É Muaig =) 
«HI 2 

7) di. = C dn 
dt ” dt 


Portanto, se se estabelecer a correspondência 
8) u<>a 


de 
9 ia | 
) E <>M, 4x 


10) <> [M, ax dt 


11) Cu Ea dp 1 
A x 


n —— 
"3 


12) Li + <> 


1) Est 


ter-se-á, para o análogo eléctrico, e com a mesma 
correspondência de grandezas, uma equação idên- 
tica a 2). 

Se, por exemplo, se quiser simular o efeito 
dinâmico de um momento M Ax aplicado na 
vizinhança da secção n da barra, entre o tempo 
t, e t,, conforme indicado na fig. 8-4a, haverá 
que, no análogo, considerar um gerador de 
corrente in , variando a corrente no tempo 
conforme indicado na fig. 8-4b (supondo igual 
escala de tempos no sistema mecânico e no 
análogo). 


M Ax 


n 


a) 


<>IM Axdt 
b) 


t, t> t 


Fig. 8-4 — Exemplo de relação entre a variação no 

tempo do momento de torção exterior, M, Ax, aplicando 

a um troço Ax da barra na vizinhança do ponto n, 

e a variação no tempo da corrente ilV correspondente, 
no análogo. 
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9. REGIMES LIVRES. FREQUÊNCIAS 
PRÓPRIAS. REPRESENTAÇÃO MODAL. 
ANÁLISE DE VIBRAÇÕES. 


Suponhamos que se pretende conhecer o regime 
livre de um sistema mecânico em determinadas 
condições, por exemplo, partindo do repouso, 
em consequência de uma perturbação impulsiva 
ou de uma perturbação escalão de determinada 
grandeza. 

Para concretizar, suponhamos o caso de uma 
viga, analisado em 7., e que se pretende conhecer 
o regime livre consequente à aplicação impulsiva 
de uma força em determinado ponto (fig. 9-1a). 
Aplicando no modelo eléctrico, análogo do sistema 
mecânico em regime dinâmico, por meió do gera- 
dor de corrente análoga dessa força, um impulso 
de corrente, têm-se reproduzidas as condições 
causadoras do regime livre. Registando, por exem- 
plo oscilograficamente, as várias grandezas eléctri- 
cas, têm-se, em função do tempo, os análogos 
das várias grandezas que caracterizam o regime 
livre. 


t t 


Fig. 9-1 — Exemplos de variação no tempo de excita- 
ções para as quais se pretendam analisar regimes 
livres: 

a) Força de tipo impulsivo, 

b) Força de tipo escalão. 


Suponhamos agora que se pretende determinar 
a resposta do sistema mecânico a uma exitação 
sinusoidal, de frequência f e pulsação w», em fun- 
ção de f, ou de 19, 

Aplicando, ao análogo eléctrico do sistema 
dinâmico, o análogo dessa excitação, tem-se, 
para determinada frequência f, por medida das 
grandezas eléctricas análogas, a resposta do sis- 
tema para esse valor de f. Variando a frequência 
podem obter-se, em função de f, as grandezas 
que se pretenda determinar. 

Tem, por vezes, particular interesse, a determi- 
nação das frequências próprias do sistema. No 
caso de sistemas não dissipativos, tem-se, para 
as frequências próprias, uma relação infinita entre 
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a «resposta» do sistema (caracterizada por deter- 
minada grandeza, por exemplo uma velocidade) 
e a «excitação» (caracterizada por determinada 
grandeza, por exemplo uma força). Como, prati- 
camente, há alguma dissipação, essa relação não 
é infinita, mas muito elevada. As frequências 
próprias corresponderão assim, muito aproxima- 
damente, a máximos relativos da «resposta», 
tendo a «excitação» amplitude constante e fre- 
quência variável. 

Seja, por exemplo, o caso representado na fig. 


9-2. As frequências próprias do sistema são F,, 
E, Eu E, 


Fig. 9-2 — Representação, em função da frequência f, 
da relação entre a excitação, caracterizada por exemplo 
por uma força com uma única componente não nula 
tendo por análogo uma corrente de amplitude 1, 
e a resposta do sistema, caracterizada por exemplo 
por uma velocidade com uma única componente não 
nula tendo por análogo uma tensão de amplitude U, 


No caso de sistemas de parâmetros distribuídos, 
a sua representação por um número finito de 
circuitos de parâmetros concentrados introduz um 
erro, equivalente à aproximação das diferenças 
finitas, erro em geral crescente com a frequência. 
Ha, assim, uma frequência máxima, para a qual, 
com determinada aproximação, é válido um deter- 
minado número finito de elementos. Em ordem 
de grandeza, há necessidade de considerar uma 
subdivisão em pelo menos cerca de 10 elementos 
por comprimento de onda, para as frequências 
do espectro das grandezas que interesse consi- 
derar com razoável aproximação. 

Nota-se que, muitas vezes, é possível aperfei- 
çoar o critério de definição e correspondência das 
zonas do sistema de parâmetros distribuídos 
representados por elementos de parâmetros con- 
centrados, reduzindo os erros resultantes de 
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quantificação, ou, o que é o mesmo, conseguindo 
igual precisão com menos elementos. 

Consideremos um sistema mecânico e uma 
parte do mesmo, constituindo um subsistema 
passivo, ligado ao resto do sistema apenas por 
um ponto de ligação e pela força exercida nesse 
ponto. 

Nessas condições, para efeito do sistema em 
conjunto, interessa apenas representar o sub- 
sistema pela relação entre a força e o desloca- 
mento ou a velocidade do ponto de ligação, ou, 
no análogo, por exemplo a relação entre a cor- 
rente e a tensão num par de terminais, na hipó- 
tese de ser suficiente considerar apenas uma 
componente do deslocamento e uma componente 
da força. 

No caso de um subsistema não dissipativo, 
para efeitos do par de terminais que, nessa hipó- 
tese, caracteriza as variáveis de ligação, o análogo 
do subsistema pode ser substituído, por exemplo, 
por um circuito do tipo do representado na 
fig. 9-3, correspondendo a capacidade C, à fre- 
quência própria «nula» (deslocamento do sub- 
sistema em conjunto sem deformação) e tendo 


Fig. 9-3 


os sucessivos paralelos indução-condensador 
frequências de ressonância correspondentes às 
frequências próprias do sistema mecânico. Desta 
forma, se for suficiente considerar n frequências 
próprias do subsistema, para efeitos do análogo 
do sistema conjunto o subsistema pode ser 
representado com um circuito constituído apenas 
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por n (ou n-1) induções e n capacidades. Essas 
induções e capacidades são definidas, a menos de 
um factor de proporcionalidade, pelos pólos e 
zeros da relação entre a resposta e a excitação 
do subsistema. Esses pólos e zeros verificam-se 
para frequências reais se o subsistema for não 
dissipativo. 

Na hipótese de não ser suficiente considerar 
apenas uma componente do deslocamento e uma 
componente da força, para a interacção entre o 
subsistema e o resto do sistema, e no caso de 
o subsistema, sendo passivo, ser actuado em 
mais de um ponto, ele pode ser considerado, 
para efeitos do sistema conjunto, por meio de um 
análogo do tipo representado na fig. 9-4, corres- 
pondendo as correntes F, às várias forças da 
interligação, a capacidade de C, à «frequência 
própria» nula, e as frequências de ressonância 
dos vários paralelos indução-capacidade às várias 
frequências próprias do sistema mecânico. Asso- 
ciam-se, directamente, desta forma, os elementos 
do análogo eléctrico aos «modos» do sistema 
mecânico. 


e — : — — .— = e — — — — Transformador | 


E: as — O — — -— — —— — Transformador ] 


Fig. 9-4 


10. ANÁLOGOS DE ELEMENTOS 
NÃO LINEARES 


Há, por vezes, necessidade de considerar rela- 
ções não lineares. Na generalidade dos casos, se 
não se pretender grande rigor na sua represen- 
tação, é possível obter análogos simples que 
traduzam a não linearidade. Para relações com- 
plexas e em que se pretenda grande rigor, a 
obtenção do análogo pode exigir soluções relativa- 
mente elaboradas ou de custo apreciável. 
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O elemento não linear mais simples e mais 
utilizado é constituído pelo díodo que, estilizada- 
mente e em determinadas condições, é carac- 
terizado por uma relação entre a corrente num 
par de terminais e a tensão entre esses terminais 
do tipo da indicada na fig. 10-1. Nalguns díodos 
são utilizáveis os ramos A e B dessa relação, 
noutros, também o ramo C. São representados 
os primeiros pela forma indicada na fig. 10-1 a, 
os segundos na fig. 10-1 b. 


Fig. 10-r — Característica ideal de um díodo, sendo i a 
corrente nos seus terminais (em sentidos opostos) e u a 
tensão entre os seus dois terminais. 


a) Representação esquemática de um díodo utilizável 
nos ramos À e B da sua característica. 

b) Representação esquemática de um díodo utilizável 
nos ramos À, Be C da sua característica. 


Apresento a título ilustrativo, três exemplos 
muito simples de análogos utilizando díodos. 


a) Mola não linear em regime estático 


Consideremos uma mola não linear em regime 
estático, com a relação entre a força F e a defor- 
mação X indicada na fig. 10-2. Em regime esta- 
cionário, fazendo corresponder a força à corrente | 
e a deformação à tensão U, ter-se-á o análogo 
indicado na fig. 10-3. 

Para tensões compreendidas entre U, e U, um 
de cada par de díodos em série do análogo 
(D eD,, e D,e D,) encontra-se no ramo B da 
sua característica (fig. 10-1), sendo nulas as 
correntes I” e 1”, e a relação entre a tensão U 
e a corrente I é definida por R,. Para tensões 
compreendidas entre U, e U, ou entre U, e U, 
verifica-se o mesmo quanto ao par de díodos 
D. e D, (I” é nula), mas não quanto ao par de 
díodos D, e D,, e a inclinação da curva que traduz 
a relação entre U e I é 
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Para corrente 1 inferior à correspondente a U, 
ou superior à correspondente a U,, a condição 


Fig. 10-2— Exemplo de relação entre a deformação, X, e 

a força F, numa mola não linear, e da relação correspon- 

dente no análogo, correspondendo a X uma tensão U ea I 
uma corrente 1. 


Fig. 10-3 — Análogo de uma mola não linear com a 
característica de fig. r0-2, em regime estático. O aná- 
logo é válido sendo U e |, análogos de X e F, os valo- 
res instantâneos da corrente e tensão (ou os valores 
constantes, em regime estacionário). 


citada deixa também de se verificar para D, e D,, 
impondo um destes díodos, segundo o ramo € 
(fig. 10-1) da sua característica, a tensão U. 
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b) Força de atrito 


Suponhamos, em regime dinâmico, uma força 
de atrito, caracterizada por uma relação entre a 
velocidade, v, e a força, f, do tipo da representada 
na fig. 10-4, 


Fig. 10-4 — Relação entre a velocidade, v, e a força f, 
do tipo de uma força de atrito, e relação correspon- 
dente entre os respectivos análogos, tensão, u, cor- 
rente i, tendo os restantes parâmetros o significado 
indicado no circuito análogo da fig. 10-5. 


Consideremos o circuito indicado na fig. 10-5, 
sendo E a tensão, estacionária, de um gerador 
de tensão. Fazendo corresponder a força f à 


x ; - E 
corrente i, a tensão u à velocidade v, = à incli- 


R 


nação da relação entre v e f, e e a f.,, tem-se 


um análogo da relação entre f e v. 


—— e 


Fig. 10-5 — Circuito análogo de uma relação velocidade- 
“força de um sistema com atrito, com a característica 
indicada na fig. 10-4. O gerador de tensão do análogo 
é um gerador de tensão É, constante no tempo, e os 
análogos de fe v são os valores instantâneos dei eu. 
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c) Relação deslocamento — deformação de um corpo 
plástico 


Consideremos o circuito representado na fig. 
10-6. 


Fig. 10-6 — Circuito tal que, em determinadas condi- 
ções e com determinada aproximação, traduz uma 
relação entre i (análogo de uma força) e F udt (análogo 
de um deslocamento) do tipo da indicada na fig. 10-7. 


Em determinadas condições, e com certa 
aproximação, este circuito é caracterizado por 
uma relação entre i (análogo de força) e fu dt 
(análogo de deformação) do tipo da indicada na 
fig. 10-7, podendo ser desta forma utilizado para 
simular deformações de tipo plástico. 


i<>F 


Ju dt<>X 


Fig. 10-7 — Relação entre uma força F (tendo por aná- 

logo o valor instantâneo de uma corrente i) e um deslo- 

camento X (tendo por análogo o valor instantâneo de 
/ udt), do tipo da de um corpo plástico. 
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11) CÁLCULO ANALÓGICO DE DEFORMA- 
ÇÃO DE PLACAS 


O cálculo de deformação estática de placas 
elásticas lineares é, em geral e se satisfeitas deter- 
minadas condições, redutível à resolução de 
equações do tipo 

2 2 2 9 
vi D+ Ly A 
dx dy” dx” 


ou semelhantes, sendo x e y as coordenadas 
cartesianas da placa, a um parâmetro, função de 
x e y (nulo nalgumas hipóteses de deformação) 
e F uma função de x e y. 


Fazendo 
») G=Guyyl= > 4 LE, 
age dy” 
tem-se 
"G |, SG mê 
dxº dy? o 
3) 2 2 
Esto 
| dx dy” 


Tem-se, assim, o cálculo da deformação de 
placas reduzido à resolução simultânea de duas 
equações de Poisson, da forma 3), o que é possí- 
vel com análogos constituídos apenas por resis- 
tências (além de fontes de tensão e-ou-corrente). 

No caso de placas sujeitas a forças e defor- 
mações segundo a direcção z, ortogonal a x e y, 
tem-se o seguinte significado dos parâmetros das 
equações 3): 


4) F==z deformação transversal 


9 
D 


5) a= 


“sendo q a densidade de carga e D um pará- 
metro de rigidez 


6) TES 
sendo M o «momento flector total» tal que 
M.+M 
» Mm OT 
D(1+p) 


sendo M, e M, momentos flectores referi- 
dos às direcções x e y, 


9“ z dz 
8) M,=-D(— +e —) 
dx” dy 
y2 3? 
9) M=—D( OZ pu SÍ) 
dy” dx 


e u a relação de Poisson. 


No caso de placas sujeitas a forças e defor- 
mação apenas segundo as direcções x e y (defor- 
mação transversal nula) tem-se 


10) a=0 


11) G=0, +69, 


9“F 
12) o Pa 
dy 
dE 
1) = 
dx 
Ni es E 
nó dx dy 


sendo 7,, Ty Txy às tensões referidas às direcções 
x, Y. 

Consideremos um quadriculado segundo as 
direcções x e y e o circuito representado na 
fig. 11-1, em que duas malhas de resistência 
R e r correspondam a esse quadriculado da 
placa, sendo a proporção das várias resistências 
e ordem de grandeza das tensões e correntes 
aplicadas e das tensões resultantes tais que, com 
suficiente aproximação : 


— as correntes i sejam desprezáveis para 
efeitos das tensões na malha de resis- 
tências R; 

— as correntes i sejam definidas pelas tensões 
da malha de resistência R e praticamente 
não afectadas pelas tensões da malha de 
resistências r; 

— as correntes | sejam proporcionais ao pará- 
metro a nos pontos x, y correspondentes 
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Fig. 11-1 — Representação esquemática de um circuito 
constituído por dois quadriculados de resistências, R er, 
correspondentes a uma decomposição num quadriculado 
segundo as direcções x e y de uma placa, sendo os nós dos 
dois quadriculados ligados por resistências ;. As tensões, 
em relação a um nó de referência, dos nós do quadriculado 
R, são o análogo dos valores da função 6, e as dos nós 
do quadriculado r são o análogo dos valores da função F, 
nas condições indicadas no texto. 


aos nós da malha de resistências R em que 
são injectadas. 


Nestas condições, dentro da aproximação de 
diferenças finitas, as tensões dos nós das malhas 
de resistências R er, em relação a um nó de 
referência, serão, respectivamente, os análogos 
de Ge F. 

Desta forma, impondo condições análogas das 
condições de carga e ligação da placa, têm-se, 
por medida directa no análogo, as grandezas 
G e F, podendo obter-se a partir delas as várias 
grandezas que interesse determinar (eventual- 
mente por um circuito «análogo» das relações 
entre G, F e essas grandezas). 
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